Realistic transport modeling for a superconducting nanowire with Majorana fermions 
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Motivated by recent experiments searching for Majorana fermions (MFs) in hybrid 
semiconducting-superconducting nanostructures and by subsequent theoretical interpretations, we 
consider the so far most realistic model (including disorder) and analyze its transport behavior nu- 
merically. In particular, we include in the model superconducting contacts used in the experiments 
to extract the current. We show that important new features emerge that are absent in simpler 
models, such as the enhanced visibility of the topological gap for increased spin-orbit interaction. 
We find oscillations of the zero bias peak as function of magnetic field and explain their origin. Even 
taking into account all the possible (known) ingredients of the experiments and exploring many pa- 
rameter regimes for MFs, we are not able to reach a satisfactory agreement with the reported data. 
Thus, a different physical origin for the observed zero-bias peak cannot be excluded. 

PACS numbers: 74.45. +c, 73.63.Nm, 74.78.Na 



The experimental search [1-3 of Majorana fermions 
(MFs) predicted to occur in condensed matter sys- 
tems [^'-n^O^ is very challenging due to the fact that these 
exotic quasiparticles are characterized by zero effective 
coupling to electromagnetic fields. Only an indirect iden- 
tification is possible via induced signatures. For instance, 
the presence of a MF is expected to induce a zero-bias 
conductance peak (ZBP) in a tunnel-spectroscopy trans- 
port experiment [TT, T?^. However, such features while 
being easily observable, are not an unambiguous demon- 
stration of MF physics. The same ZBPs can indeed be 
induced by various different mechanisms, including the 
Kondo effect [13 , Andreev bound states [I4, weak an- 
tilocalization and reflectionless tunneling [15] . Ruling out 
these alternative scenarios has been made difficult by the 
fact that the measured conductance has a rather unex- 
pected behavior, so that also the explanation in terms of 
MFs is not entirely satisfactory. 

A typical experimental setup [1-3 (see Fig. [T]) con- 
sists of a spin-orbit-coupled semiconducting nanowire de- 
posited on or coated with a bulk 5-wave superconductor 
on one end and contacted through a tunnel barrier by a 
normal lead, on the other end. Part of the nanowire it- 
self is assumed to be in a superconducting state, induced 
through proximity effect by the bulk superconducting 
lead. 

The magnetic-field-induced transition to the topolog- 
ical phase in the nanowire is accompanied by a closing 
and reopening of the excitation gap [5HlQ]. The topo- 
logical phase formally persists for all values of magnetic 
field above a critical in a one-band model, while it 
could have a finite upper critical field in a multi-band 
model, where bands are crossing at large Zeeman split- 
tings and hybridization of MFs can take place. How- 
ever, in experiments one typically explores magnetic field 
regimes where only one band undergoes a topological 
transition [iHS]. Since the bulk superconductor always 
remains in a non-topological state, the topological sec- 
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FIG. 1. The schematics of the NSS^ geometry setup we con- 
sider in this work- (top panel). The sample under investiga- 
tion is connected on the left to a semi-infinite normal lead (N, 
orange) and on the right to a semi-infinite bulk s-wave super- 
conducting lead (S, green). It consists of a normal nanowire 
section (NW, gray), where a potential barrier U{x) (black) is 
created, and a proximity- induced superconducting nanowire 
section (SW, gray). We allow for static disorder w{x,y) (red 
crosses) in the nanowire. The spatial dependence of all the pa- 
rameters entering the Hamiltonian in Eq. ([T]) is qualitatively 
depicted in the bottom panel. 



tion has a finite length L^, and MFs localized at its ends 
can overlap, depending on their localization length ^m- 
Given that the latter depends also on the value of mag- 
netic field, one expects to observe a ZBP that splits for 
large enough 5-fields. 

The experiments conducted so far [T]-[3] reported fea- 
tures which are partially consistent with the existence 
of Majorana end-states in the nanowire, namely a (non- 
quantized) ZBP when a magnetic field of sufficient 
strength is applied along the nanowire. However, quan- 
titative agreement with the theory is still missing, and in 
particular the following points have to be clarified: 

a) The most evident discrepancy between experiment 
and theory is the absence of any experimental signa- 
ture related to the closing of the excitation gap in the 
nanowire. Recently, this fact has been ascribed to the 
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spatial distribution of the wave-functions associated with 
the closing-gap-states in the regime of low chemical po- 
tential [lillT]. 

b) The ZBP in the experiments appears above a certain 
magnetic field, persists over a finite range of values, and 
then disappears again symmetrically, rather than mutat- 
ing into a pair of split peaks as expected for MFs. 

c) The ZBP conductance is not quantized, with values 
being a small fraction of the quantum of conductance 
2e^ /h [l]-[3], while in theory the Majorana is providing a 
full transport channel and produces a quantized conduc- 
tance peak G = 2eV/i dHIIlITl]. 

d) The proximity-induced gap in the nanowire de- 
pends only very weakly on magnetic field in the mea- 
sured dl/dV curves, and the corresponding conductance 
signal decreases significantly above a certain magnetic 
field [U |3] . Such sudden reduction is not predicted by the 
theory, and the gap should close much faster as a func- 
tion of magnetic field. This issue has not been pointed 
out in previous theoretical studies. 

In an attempt to find explanations for above issues, 
we have performed numerical calculations of the two- 
terminal conductance in an hybrid structure with a su- 
perconducting spin-orbit-coupled nanowire contacted on 
one side by a normal metallic lead and on the other side 
by a standard s-wave superconducting lead, referred to 
as NSS' setup (see Fig. [T]), which closely resembles the 
actual geometry of the experiment. The conductance is 
calculated within the standard framework of scattering 
theory [20], with the help of the recursive Green's func- 
tion technique ^ applied to a real-space tight-binding 
Hamiltonian. Such a procedure allows us to model a 
complex structure like the one effectively involved in an 
experimental situation that is not amenable to analytical 
approaches [22 . 

To be specific, we focus now on energy scales appro- 
priate for InSb nanowires [H , with the exceptions and 
specifications described below. Nonidealities and compli- 
cations such as multiple occupied sub-bands, presence of 
disorder, finite width of electrostatic barriers, finite co- 
herence lengths of superconducting pairing amplitudes, 
and non-zero temperature are taken into account. 

Our study reveals important features not emphasized 
so far. In this respect, the presence of the bulk supercon- 
ductor turns out to be decisive. We briefiy summarize 
here our main findings. In some regimes, the closing of 
the gap becomes visible in the conductance signal, while 
it does not in an NS-setup simulation. Further, the gap- 
edge conductance peak value decreases as a function of 
magnetic field, a feature which is also not captured by 
simpler models. We find oscillations of the zero bias peak 
as function of B and explain their origin. We argue that 
disorder is unlikely to be the explanation for the ZBPs 
observed in experiments. Further, we show that the tun- 
nel barrier plays an important role in determining which 
features are actually visible. Finally, according to our 



results the experimental dl/dV behavior seems to point 
to a SOI strength larger than the one reported. 
Model. We consider a two-dimensional rectangular 
nanowire of length L along the ;r-direction and lateral ex- 
tension W in the y-direction. The general tight-binding 
Hamiltonian describing the different sections of the setup 
in the presence of a magnetic field has the form 



m 

E^- (4,t4,i + h.c.), 



Cm, (3 
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where t is the hopping amplitude and a is the spin-fiip 
hopping amplitude, related to the spin-orbit interaction 
(SOI) energy E^o through the relation E^o = o? jt. The 
sum in the first row runs over all lattice sites m and 
relative nearest neighbors (m + d). Implicit summation 
over repeated spin indices is assumed. The constant /io 
is chosen to set the common chemical potential to the 
zero-field bottom of the topmost band in the nanowire, 
and its value depends on the number of nanowire sub- 
bands (i.e. on W). The term e^n = — /im + + 
accounts for local variations of the chemical potential, for 
the tunnel-barrier potential /7m, and includes an on-site 
random potential Wm which models Anderson disorder. 
The tunnel barrier is assumed to have a Gaussian pro- 
file with height and width A. An external magnetic 
field B is oriented along the nanowire axis {x) and in- 
duces a Zeeman splitting 2V^ = ^m^B^- Finally, A 
is the pairing amplitude and can either account for the 
native superconductivity in the bulk 5-wave supercon- 
ducting lead (Aq) or for the proximity-induced pairing 
in the nanowire (A^). All these quantities are taken to be 
site-dependent along the x direction (except Wm^ which 
is taken to be completely random), so that by varying 
their value we can model physically different parts of the 
setup. The normal lead is characterized by a = 0, /i :^ /io 
(i.e. metallic regime), ^ = 2, Wm = 0, and A = 0. The 
nanowire is characterized by finite a = aR, chemical po- 
tential /i close to the bottom of the top-most band, 
^ = 50 appropriate for InSb nanowires, and A varying 
from in the normal section to A^ in the proximized 
section. The nanowire is adiabatically connected to a 
metallic superconducting lead with a = 0,/i:^/io,^ = 2, 
= 0, and A = Ao > A^. By adopting the sequential 
geometry of Fig. [l] we slightly simplify here the actual 
situation in the experiments [IHSl? where the nanowire 
is side-contacted, or top-contacted, and the current does 
not follow a straight path. 

In a simpler model the nanowire can be taken to 
extend to infinity, without any external superconductor, 
in an effective NS geometry. In our formalism this 
corresponds to taking the superconducting lead to be 
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physically identical to the nanowire, with a single pairing 
amplitude A^. In such a configuration, the second Majo- 
rana end-state is always moved to infinity, and the ZBP 
is locked to zero for all values of magnetic field B > B^^ 
whereby the topological transition occurs at the critical 
field {g/j.B/2)Bc = ^A^H^ [SHIO]. We wih sometimes 
switch to this NS configuration in order to make contact 
to previous theoretical studies [161 [TTl [23H27] and to 
understand the effect of the bulk superconductor. 




FIG. 2. Effect of larger SOI strength. We plot here the dif- 
ferential conductance dl/dV evaluated as a function of bias 
voltage V and Zeeman energy Vz . Panels a) and b) refer to the 
NS configuration, while c) and d) refer to the NSS^ setup. The 
parameters used here are = 250 /xeV, Ao = 2.1 meV (only 
NSS^), Uo — 45 meV, A = 1 nm (narrow barrier), Ln — 0, 
= 6 /im (only NSS^) and /x = 0, which corresponds to a 
critical = A^. For the case of InSb, the plotted range 
Vz = — 6A^ corresponds to 5 = — 1 T. Temperature is 
set to T = 75 mK. (Left column) Eso — 50 /xeV. (Right col- 
umn) Eso = 200 /ieV. Larger SOI yields a slower closing of 
the A;F-gap Afcp(B), in both configurations, where /cp is the 
Fermi momentum. Notice that in the NSS^ case the A;F-gap 
signal decreases in intensity as the magnetic field grows. 

Let us first note that the actual value of the SOI a in 
the experiments is not known, as also noticed in Ref. [28] . 
since the only available measurements have been per- 
formed in a different setup, where the effective Rashba 
SOI could be largely modified. Similarly, the value of 
the proximity-induced pairing amplitude is not directly 
accessible, and one can only infer its value from the con- 
ductance peaks. Having said that, it becomes interest- 
ing and even necessary to consider regimes with different 
SOI, or different proximity-induced gap magnitudes. 

The first important point we want to make is that by 
assuming that the actual SOI is larger than the reported 
one (e.g. E^o = 50 /ieV in Ref. [T|), one can get a sub- 



stantial improvement in the calculated dl/dV behavior, 
with features resembling closer the ones observed in ex- 
periments [IHS]. In other words, the data gives hints of a 
stronger SOI. In particular, we note the following facts. 

1) Under the assumption that the ZBP observed in 
Refs. [IHS] arises from MFs, one must conclude that 
the chemical potential in the topological region is low, 
jj. 0^ since the ZBP emerges at low magnetic field, 
^q/JjbB A^ for g = 50, which implies a low fi in view 
of the relation between Bc^ /i, and A^. However, such a 
low /i, together with the reported SOI values pL^, would 
generate a rapid closing of the /cp-gap Af^^ as a func- 
tion of magnetic field. This is indeed what we observe in 
our numerical calculations when we work in the regime 
/i 0, Eso = 50 /ieV, both in the NS and NSS' setup, see 
Fig. |2|l) and c) , respectively. Note that in the NS case 
the ZBP stays at zero bias for all magnetic fields, whereas 
in the NSS' case the finite length of the topological re- 
gion induces an oscillating splitting, analyzed below. In 
the same figure we show that a stronger SOI gives better 
agreement with the measured A/ep(5), both in the NS 
setup j^, see panel b), and in the NSS' setup, shown in 
panel d). Note that this latter SOI effect is independent 
of the nature of the observed ZBP. 
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FIG. 3. Effect of larger SOI strength on disorder, NSS' case. 
The parameter values are the same as in Fig. |2] In addition, 
a realistic disorder Wm G [—2.4, 2.4] meV (corresponding to 
a mean free-path imfp — 2.5 /xm [28 ) is included over the 
entire nanowire length L ~ 2.5 /im. We do not average over 
disorder configurations, a) Eso = 50 /xeV. b) Eso = 200 /ieV. 
In the weak SOI regime, the disorder lowers or destroys the 
gap from other sub-bands, bringing many supra-gap states 
down close to the Fermi level, where they cluster in some cases 
into a finite-extension ZBP, like in panel a). Such clustering 
is, however, removed for stronger SOI [28], see panel b). 

2) When realistic Anderson disorder is included in the 
model, the closing of the gap becomes visible again even 
in the jii regime [T71 |26l |28] , re- introducing a discrep- 
ancy with the experimental observations [1^3]. Disorder 
in a weakly spin-orbit-coupled nanowire causes a number 
of sub-gap states to appear, some of which cluster around 
zero-energy and possibly give rise to a non-topological 
ZBP, more markedly in the case of finite ja [28]. Such 
states are coming from other sub-bands, for which the 
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effective mini-gap gets reduced in the presence of disor- 
der. This is substantiated by the fact that the ZBP in 
Fig.|3^) has a conductance peak value larger than 2e^//i, 
implying that it cannot come from the Andreev signal of 
a single state. In the case of stronger SOI, the effect of 
disorder is partially suppressed, and less sub-gap states 
(or no states at all) are observed (see Fig. [3| , being more 
compatible with experiments pii3 . Thus, disorder is un- 
likely to explain the observed ZBPs. 
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FIG. 5. Schematic dependence of the MF localization length 
on magnetic field B. According to the theory for a one- 
band semi-infinite nanowire 9l , a MF emerges when the 
magnetic field exceeds a critical value Be = Al + /j^'^/q/jb, 
and the system goes from the non-topological (grey) to the 
topological (yellow) regime. However, for a nanowire of finite 
length L^, due to overlap of the MFs from each end, the 
additional condition for the observation of a MF is < 
. Considering typical dependences of on magnetic 
field [29 , we predict that the MF should be observed for B* < 
B < Bi\ where ^yiiBl) ^ ^uiBD ^ L./2 (cf. Fig. [If. 



FIG. 4. Same parameters as in previous figures, without dis- 
order and at T = 0. (a) = 50 /xeV. Note the oscillations of 
the ZBP, for explanations see text, (b) = 200 /xeV. With 
large enough SOI strength the gap-closing becomes visible, 
even for the considered case of /i 0, while it is not visible 
in an NS setup with the same parameters. The dl/dV peaks 
coming from the gap-closing have, however, very small width 
and they get washed out by realistic temperatures. 

3) As a consequence of finite-length of the topologi- 
cal section (L^) and of the 5-dependence of /cp, we ob- 
serve that the ZBP splitting exhibits oscillations as B is 
swept, see Fig. [4^). To explain this, we recall that in 
the weak- SOI limit the MF wave function has an expo- 
nentially decaying envelope with localization length 
and a fast-oscillating part ~ sm{kYx) [29 . If the mag- 
netic field exceeds a critical value B** (see Fig. [5|, the 
two end-MFs overlap and split away from zero energy. 
Since increases with B [5 -10, 29 , so does the split- 
ting. However, if k^L^ becomes an integer multiple of 
TT as function of 5, the ZBP splitting returns to zero, 
leading to oscillations with period given by 

where m is the band mass and a the lattice constant. Us- 
ing parameter values corresponding to Fig.jlj t/A^ = 40, 
L^/a = 200, we obtain quantitative agreement with the 
simulated ZBP oscillations. Since the critical field 5** 
increases with SOI [29], the ZBP splitting and related 
oscillations occur at larger fields, see Fig.[4]3). Note that 
these oscillations are quite robust against temperature 
effects, see Fig. This behavior of the ZBP is quite 
remarkable and provides an additional signature to iden- 
tify MFs experimentally. 



We note in passing that in the NSS' setup the SOI 
affects the visibility of the gap-closing, see Fig. [4j This 
is not the case for the NS setup. 

Next we address further issues that have received less 
attention in the literature so far. 

4) The position of the proximity gap Akp{B = 0) de- 
duced from the dl/dV curves is in general different from 
the nominal value inserted in the microscopic Hamil- 
tonian, Eq. ([T]). This observation is important, since it 
could mean that deducing the value of the proximity- 
induced gap from the conductance spectra is not a cor- 
rect procedure [III^. Such an energy shift can be due to 
the presence of a finite normal region between the tunnel 
barrier and the effective NS interface. In that case, the 
observed peak occurs at a bias voltage lower than the 
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FIG. 6. Role of the tunnel barrier U{x). The parameters are 
chosen as in the previous figures. Here we kept T = in 
order to show the effect of the barrier smoothness alone, a) 
Gaussian tunnel barrier with width A = 1 nm. The Majorana- 
induced ZBP is fully visible, with maximal weight dl/dV — 
2e^ /h at the largest magnetic fields. The closing of the gap 
is, however, nearly absent, b) Same system, but with A = 50 
nm, a value closer to the experimental situation. Gap-closing 
and ZBP are completely absent. 
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nominal A^. By increasing the distance of the barrier 
from the NS interface, one can move the conductance 
peak deeper inside the gap and eventually even intro- 
duce additional peaks, similar to the McMillan- Ro well 
resonances [30^. Alternatively, in the NSS' configuration 
the peak corresponding to itself can be viewed as a 
sub-gap resonance of the larger gap Aq, and its position 
can be changed by varying the distance of the tunnel bar- 
rier from the S-S' interface. In this case, the peak has 
energy larger than the microscopic value of A^. 

5) In both NS and NSS' configurations, the tunnel bar- 
rier plays an important role — it determines the transmis- 
sion of each transport channel, which in turn sets the 
width of the McMillan- Rowell resonances [30] . When the 
resonance width becomes smaller than the temperature, 
the resonance is practically invisible. Consequently, if 
the barrier is wide enough, no sub-gap features at all are 
present in the dl/dV signal. This can be attributed to 
the associated maximum momentum mismatch ~ 
that the barrier can provide to an incoming electron. If 
the tunnel barrier is chosen to be sharp, all the states 
present in the nanowire could in principle be made visi- 
ble. However, that is not a realistic choice, since a typ- 
ical barrier in experiments has a characteristic width of 
~ 50 nm. For such values we already observe a momen- 
tum filtering due to the barrier, leading in some cases to 
a complete disappearance of MF signatures, see Fig. |6] 
Again, introducing disorder can make the aforementioned 
sub-gap features re-appear. Therefore, it is the combined 
effect of barrier shape, SOI, and disorder strength that 
determines the visibility of MFs. 

In conclusion, by numerically simulating a more com- 
plex and realistic setup than before, we have obtained 
new features in the transport simulations, which are sim- 
ilar to the ones observed in experiments. However, even 
after considerable effort, we did not manage to repro- 
duce all such features in a single configuration, and we 
still lack a satisfactory agreement with experiments. In 
particular, the exact shape of the measured ZBP is not 
very compatible with the picture of MF states that form 
and then split as a function of magnetic field. Thus, a 
different physical origin for the observed ZBP [T]-[3] can- 
not be excluded. More precisely, from our findings it 
seems possible that in the experiments the MF-features 
are essentially invisible and the observed ZBP is coming 
from some different co-existing phenomenon, like Kondo 
effect, which seems indeed to yield a similar behavior [31 . 
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